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An interplay between charge discreteness, coherent scattering and Coulomb interaction yields 
nontrivial effects in quantum transport. We derive a real time effective action and an equivalent 
quantum Langevin equation for an arbitrary coherent scatterer and evaluate its current-voltage 
characteristics in the presence of interactions. Within our model, at large conductances Go and low 
T (but outside the instanton-dominated regime) the interaction correction to Go saturates and causes 
conductance suppression by a universal factor which depends only on the type of the conductor. 



Coulomb effects in mesoscopic tunnel junctions have 
recently received a great deal of attention jl]-|]]. One of 
the remarkable features of such systems is that charge 
quantization (and, hence, Coulomb blockade) persists 
even for junctions with low resistances 1/Gf <C Rq = 
h/e 2 w 25.8 kil. In this limit an effective Coulomb gap 
Eq for a junction with the "bare" charging energy Eq 
suffers exponential renormalization [^) 



Ec/Ec <x exp(-G t i? Q /2), 



(1) 



but remains finite even at very large values of GtRQ- 
Eq. (0) was confirmed in several later studies both ana- 
lytically and numerically |@,|J. Experiments clearly 
demonstrated the existence of charging effects for the val- 
ues of Gt as large as GtRQ « 32 |J. 

Recently another interesting prediction was made by 
Nazarov [flPf] , who argued that features of charge quan- 
tization may also persist in arbitrary conductors includ- 
ing, e.g., disordered metallic wires with g = GqRq ^> f . 
Here and below G = 1/ R = (2e 2 /h)J2 n T n is the con- 
ductance of an arbitrary scatterer and T n are the trans- 
missions of its conducting modes. If one accounts for the 
spin degeneracy, the renormalized Coulomb energy for a 
general conductor derived in ]10|| takes the form 



E c /E c ^Y[R ri 



(2) 



(Eq = 0) even at T = if at least one of the conducting 
channels is fully transparent R n — jL2j. 

In this Letter we will analyze an interplay between 
Coulomb effects and quantum transport at energies 
larger than Eq (0)- We will derive a real time effective 
action and formulate a quantum Langevin equation for 
an arbitrary (albeit relatively short) conductor. At tem- 
peratures or voltages above Eq we will obtain a complete 
I — V curve at large enough g. We will demonstrate that 
Coulomb interaction leads to (partial) conductance sup- 
pression with respect to its "noninteracting" value Go- 
This suppression effect is controlled by the parameter 



(3) 



well known in the theory of shot noise |l3[ . The parame- 
ter (3 (||) equals to one for tunnel junctions and to 1/3 for 
diffusive conductors. In contrast to Eq (§), it vanishes 
only if all the conducting channels are fully transparent. 

We identify three different regimes for the interaction 
correction to Go- Let us display the results for a linear 
conductance G{T). At T/Ec 3> max(l,g) perturbation 
theory in Eq (or in 1/T) is sufficient. It yields 



Ec _/3C(3) 1 \ (Ec 
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where R n = 1 — T n . In particular, for diffusive con- 
ductors, similarly to eq. (|l|), one finds (!(]] Ec/Ec oc 
exp(— ir 2 g/8). The same result (||) follows from the effec- 
tive action derived in ^,[nj for metallic contacts within 
the quasiclassical Green functions technique. Hence, one 
can expect the effective actions PJTl|] and JToJ to be 
equivalent, perhaps up to some unimportant details. 

Eq. (j^) sets an important energy scale for the prob- 
lem in question: at temperatures below an exponen- 
tially small value Ec a conductor with g ^> 1 should 
show insulating behavior due to Coulomb effects. On the 
other hand, at larger temperatures/ voltages this insulat- 
ing behavior should not be pronounced. Furthermore, ac- 
cording to (0) Coulomb blockade is destroyed completely 



Here £(3) ~ 1.202 and g needs not to be necessarily large. 
For j > 1 there exist two further nonperturbative in 
the interaction regimes. At intermediate temperatures 
gE c exp(-g/2) «T« gE c we have 



G_ 

Go 
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where 7 ~ 0.577. Here energy relaxation plays an im- 
portant role turning the power law dependence (Q) into a 
much slower one (^) . Finally, at even lower temperatures 
T < gEcexp(—g/2) (but T > E c ) relaxation processes 
yield complete saturation of G(T): 



G/G a a\-P + 0<fi/g). 



(6) 
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It is remarkable that the result (Jfy does not depend on 
the charging energy Ec at all. In the tunneling limit (all 
T n <C 1) the regime @ does not exist. Two other regimes 
are already known for tunnel junctions: by setting /3 = 1 
in eqs. (Q), (|5|) we recover the results [^4|-p6[. 

The model and effective action. Now let us proceed 
with the derivation of the above results and the I — V 
curve. Our framework is very similar to that of Ref. Jio| . 
We will consider an arbitrary scatterer between two big 
reservoirs. The scatterer length is assumed to be shorter 
than dephasing and inelastic relaxation lengths, so that 
phase and energy relaxation may occur only in the reser- 
voirs and not during scattering. Coulomb effects in the 
scatterer region are described by an effective capacitance 
C. The charging energy Ec = e 2 /2C, temperature T 
as well as other energy scales are assumed to be smaller 
than the typical inverse scattering time (e.g. the Thou- 
less energy in the case of diffusive conductors). 

Quantum dynamics of our system is fully described 
by the evolution operator on the Keldysh contour. The 
kernel of this operator J may be represented as a path 
integral over the fermionic fields. Performing a standard 
Hubbard-Stratonovich decoupling of the interacting term 
in the Hamiltonian enables one to integrate out fermions. 
Then the kernel J acquires the form of the path integral 
over the Hubbard-Stratonovich fields on the forward (Vi) 
and backward (V 2 ) parts of the Keldysh contour 



J = J VV 1 VV 2 exp(iS[V}), 



(7) 



where S[V] is the effective action defined as 



iS[V] = 2TrlnGy 1 +i— J dt'{Vj 



LRl 



Vl m ), (8) 



where Vlr% = Vu — Vm are the voltage drops between 
the reservoirs. The Green-Keldysh matrix Gv(X\, X 2 ) 
(here X — (t, r)) obeys the 2x2 matrix equation 



iAl-£ ( ri )i + eV(Xi) 
at\ 



Gy = 5(Xi - X 2 )cr z , 



(9) 



and similarly for G12 and G 22 . Here and below integra- 
tion over the spatial coordinates is implied in the prod- 
ucts of operators. In (noh we have defined 



Uj(ti,t 2 ) — Texp 



*2 

-i J dt' (H -eVj(t',r) 

ti 



(11) 



as the evolution operators jl7| and p(t) as the density 
matrix. The latter satisfies an exact equation |l7j 

i^p = [H ,p] - (1 - p)eV 1 p + peV 2 (l - p), (12) 

with the initial condition p(t = 0) = p 0l where po is the 
equilibrium density matrix for noninteracting electrons. 

Next we define the conducting channels in a standard 
manner. They are just the transverse quantization modes 
in the reservoirs. Describing the longitudinal motion 
within one channel quasiclassically we define the free elec- 
tron Hamiltonian in the reservoirs as follows 



d_ 



(13) 



where to, n are the channel indices and v m is the chan- 
nel velocity. In every channel the coordinate y runs from 
— 00 to for the incoming waves, and from to +00 
for the outgoing ones. The scattering matrix S, which 
is assumed here to be energy independent, relates the 
amplitudes of incoming and outgoing modes as follows 



tpm(y = +0) = S„ 



(14) 



where Hq(t) is a free electron Hamiltonian for the system 
"scatterer + reservoirs", V is the diagonal 2x2 matrix 
with the elements Vy = and <r z is the Pauli matrix. 
In the last term of eq. (||) we already made use of our 
model and assumed that the fields do not depend 
on the coordinates inside the reservoirs, i.e. for the left 
(right) reservoir we put Vj(t',r) = VL(R)j(t'). 

The elements of the Green-Keldysh matrix Gy can be 
expressed as follows: 

Gn (ti M) = -h)U x (*! ,t 2 ) + iU x (*! , t)p{t)U x (t, t 2 ), 

G 21 (h,t 2 ) = -iU 2 (h,t)(l-p{t))Ui(t,t 2 ), (10) 



where S mn are the elements of the scattering matrix 
S defined in the basis ^o,m = e lky I 'y/vm- The factor 
\J v n /v m appears in ( |l4|) since we work in the basis of 
the eigenf unctions of (13) ^o,m = e lky . Finally, the 
matrix elements of the fluctuating voltages Vj(t) are: 
Vj y mn{t) = Vj : m(t)S mn , where Vj, m (t) = Vtj(i) for the 
left channels and Vj iTn (t) = Vftj(t) for the right ones. 

With the aid of Q, @ the evolution operators Q 
can be evaluated exactly. In this paragraph we will sup- 
press the Keldysh index for simplicity. Denoting the wave 
function at some initial time t\ as tl> n {t\, y), at some other 
time t 2 we find 

1pn(h,y) = i>n(tl,y- V n (t 2 - *l))Xnn(*2,*l)) (15) 

for y < or y > v n (t 2 — t\) and 



ipn{h,y) = ^ S n k* — tpk [h, —y ~ v k {t 2 - tf) 



x Xnk{t 2 M)xnk{t 2 - y/v n ,t 2 - y/v n ), 



(16) 
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for < y < v„(t 2 - h). Here we defined Xnk{t 2 ,h) = 
ex.p(i(p„(t 2 ) - ipk(h)) and ip n (ti) = f Q % dt' eV n (t'). 
On the other hand, by definition we have 

V>n(*2,Z/2) = ^ dyiU nk (t 2 ,t 1 ;y 2 ,yi)-)jjk(t 1 ,y 1 ). (17) 



Comparing ( Jl7| ) with ([l5|), (i6) and introducing a new 
coordinate r = y/v n we obtain 



U(t 2 t 1 \T 2 T 1 ) = 6{t 2 -ti -t 2 +h)e 



iip(t 2 ) 



+ 9(T 2 )e(- n )e 



-iip(t 2 ~T 2 ) 



[5-1} 



-ic^(ti) 



. (18) 



where e 1 ^ is the diagonal matrix with the elements e %Vn . 
Restoring the Keldysh index in ([l8]) we arrive at the de- 
sired result for the evolution operators U\^ 2 . 

In order to evaluate the density matrix p in the pres- 
ence of interactions one should solve a nonlinear equation 
( |l2| ) for arbitrary realizations of the fluctuating fields V\ 
and V 2 . In general this task cannot easily be accom- 
plished. Fortunately it suffices for our present purposes 
to find the density matrix for the case V\ (t) = V 2 (t) only. 
In this case eq. ([12h is trivially solved and we get 



p(t) = U(t,0)p o U(0,t), 



(19) 



where U is defined in eq. (|18[). 

In order to proceed we will make use of the quantum 
Langevin equation approach . In the case of metal- 
lic tunnel junction this approach was developed in Refs. 
0,|O|J14|. Let us define ecp+(t) = (V Lm (t) + V LR2 (t))/2 
a,nd eip~(t) = VLm{t) — VL R2 {t). The key step is to treat 
quantum dynamics of the 1^-fields within the quasiclas- 
sical approximation, i.e. to assume that fluctuations of 
p~ (t) are sufficiently small at all times. This assumption 
allows to expand the exact effective action in powers of 
p~ while keeping the full nonlinear dependence on the 
"center-of-mass" field p + . This approximation is known 
fjOjfLij to be particularly useful in the limit g ^> 1. 

Expanding TrlnGy 1 up to the second order in p~ we 
obtain 



2Tr In G 



v 



2TrhiG~ 



\ v -=o 



iS R — Si, 



(20) 



where iS R = Tr{(G n + G 22 )ip }, Si = 
Tr{Gi 2 ip G 2 \<p } and ip is the diagonal matrix with 
the elements <p>~. The zero order term in the expansion 
( ^p| ) vanishes. Evaluating the first order term iS R , with 
the aid of eqs. (10), ( |l8| ) and ( ^9|) one finds at sufficiently 
long times t 



As before, g — 2tr[£ + f] is the dimensionless conductance 
of the scatterer and t is the transmission matrix. An 
analogous calculation of the second order term Si yields 



q _ 9 

47T 2 



t t 

dt' J dt"a{t' -t")<p-(t')tp-(t") 
o o 



x{Pcos[p + (t')-p + (t")} + l-f3}, 



(22) 



where we defined a(t) = (ttT) 2 / sinh 2 [7rT£] and fig 
2tr[t+f(l - i+i)]. Combining the results (|l]) and (|2^ 
with the last term of eq. (||) we arrive at the final ex- 
pression for the effective action 



iS = i j di 
o 



C . , . _ I x 
— tp^tp H ip 



iS R -Si. (23) 



In (|2^) we also included the term which accounts for an 
external current bias I x . 

Quantum Langevin equation. The action (pl|)-(p3|) has 
the same form as one for a tunnel junction with the con- 
ductance (3g shunted by an Ohmic conductor (I — (3)g. 
Eqs. (|2l|)-(|2^) are equivalent to the Langevin equation 



C 



-ip- 



eR 



Ix = £iCOS^+ +£ 2 sin ip + + £ 3 , (24) 



where the terms in the right-hand side account for the 
current noise and are defined by the correlators 

&(*&(0)> = ~^{t) (V + 5 j2 + ^*ia) . (25) 



In the small transparency limit eqs. (|24|), (25) reduce to 
those derived before for metallic tunnel junctions fll~9| , p0"| . 
If we decompose p + (t) = eVt + 5p + (V is the average 
voltage across the conductor) and neglect the fluctuating 
part of the phase 8p + we will immediately reproduce the 
well known results |l5| ] for the current noise in mesoscopic 
conductors. 

I- V curve. In order to study the influence of Coulomb 
effects on the current- voltage characteristics for an arbi- 
trary scatterer we will make use of the exact identity 



m + m 3S\ip + , ip 

5p (t) 



JS[<p + ,<p ] = Q 



(26) 



Evaluating this path integral we set cos[ip + (t') — p + (t")] = 
cos[eV (t' — t" )] in the exponent of ( p6| ) but retain the full 
nonlincarity in SS/Sp~ . This approximation works well 
provided either g ^> 1 or max(T, eV) ^> Ec- A straight- 
forward calculation then yields 



iS R 



i_9_ 

2tt 



dt' p-(t')p+{t'). 



(21) 



X ~ R 7T 



dta(t)e- Fit \l 



')sm[eVi\, (27) 
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i r +ca 

Fit) = — / dt'a(t') {j3cos[eVf\ +1-0) 

9 J-oo 



\t' -t\ - \t'\+RC 



\/RC 



\t'\/RC 



(28) 



Eqs. (p7|) , ( p8| ) represent the central result of this paper. 

Single scatterer. In the limit g > 1 and max(eV r , T) ^> 
gEc exp(— g/2) the integral in ( |27| ) converges at times 
for which -Fft) is still small and can be neglected. In this 
limit eq. (E7h yields 



Ir 



V 
R 



— e(3Tlm 



*( 1+ f) 



- iv^ 1 



IV 

~2 



(29) 



where ^(x) is the digamma function, w = u + iv, 
u = gE c /TT 2 T and v = eV/irT. At T -> from @ 
we obtain 



fl_ = l--In|l- 



1 



.g \ (eViZC) 2 , 
while in the limit eV/Ec ^> max(l,(jf) we find 



RI* 



V-/3e/2C. 



(30) 



(31) 



For (3=1 the results (|30|), ( |3l| ) reduce to those derived 
in Refs. [ p0|pl] ]) for tunnel junctions. Eq. ( |3l| ) demon- 
strates that at large V the I — V curve of any relatively 
short conductor should be offset by the value (ie/2C due 
to Coulomb effects. For instance, in disordered conduc- 
tors this offset is expected to be only 3 times smaller than 
for a tunnel junction with the same Eq- 
At V -> from (p9h we get 



G/Gq = 1 



2£ 
9 



7 + *(l+2) + 2*( 1+ 2)]' (32) 



which yields eqs. (^) and (||) in the correspond ing limits. 
[The term with 1/15 in (0) is recovered from (p7|), (Eq).] 
In the limit max(eV, T) < gEc exp(—g/2) the integral 
(^7]) converges at very long times and the function F ( ps| ) 
cannot be disregarded. Evaluating ( p8|) at t ^> 1/RC 
we find F(t) ~ {2/g)(ln(t/RC) + 7) and performing the 
integral in ( p7|) for j > 1 we arrive at the result (|^) 
G/Gq = X/n^n- Hence, at extremely low T the inter- 
action correction to Go saturates due to Coulomb and 
relaxation effects. For diffusive conductors eq. @ yields 
G/G a 2/3. 

In the limit g 3> 1 our results are valid except for 
exponentially low T,eV < in which case instan- 

ton effects [p||^JTo|] gain importance and eventually turn 
a conductor into an insulator at T = 0. These effects 
are beyond the scope of the present paper. For tun- 
nel junctions the regime ^ cannot be realized since in 
that case Eq/Ec oc exp(— g/2). In other cases, however, 
Ec/Ec <C exp(— g/2) and the saturation of Gq(T) be- 
comes possible. Furthermore, if the instanton effects are 
suppressed (Eq — * 0), our results should apply down to 
zero temperature and voltage. 



Two scatterers. The effects discussed here can be con- 
veniently measured e.g. in the "SET transistor" configu- 
ration HH of two scatterers (such as, e.g., quantum point 
contacts) connected by a small metallic island. With 
simple modifications our results hold for such two scat- 
terer systems as well. For instance, Go is defined by eq. 
( |32| ) where R is now a sum of two resistances Ri + R2, 
« -» (.91 + 92)Ec/tt 2 T and 

P _^ ffiffg + fogi 
9 9i+92 

The I —V curve is offset at high voltages as in eq. ( |3l| ) 
with (3 — > /3i + /?2 and G being the total capacitance of 
the device. Gate modulation effects can also be treated 
along the same lines as it was done in Ref. |H . 

In summary, we studied the effect of Coulomb interac- 
tion on the I — V curve of a coherent scatterer. At low T 
its conductance is suppressed by the universal factor (||). 

We would like to thank Yu.V. Nazarov and G. Schon 
for useful discussions. 
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